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Abstract 

We present the results of spectroscopy of Lyman Break Galaxies (LBGs) at z ~ 5 in 
the J0053+1234 field with the Faint Object Camera and Spectrograph on the Subaru 
telescope. Among 5 bright candidates with z' < 25.0 mag, 2 objects are confirmed 
to be at 2 ~ 5 from their Lya emission and the continuum depression shortward of 
Lya. The EWs of Lya emission of the 2 LBGs are not so strong to be detected 
as Lya emitters, and one of them shows strong low-ionized interstellar (LIS) metal 
absorption lines. Two faint objects with z' > 25.0 mag are also confirmed to be at 
2; ~ 5, and their spectra show strong Lya emission in contrast to the bright ones. 
These results suggest a deficiency of strong Lya emission in bright LBGs at z ^ 5, 
which has been discussed in our previous paper. Combined with our previous spectra 
of LBGs at 2; ~ 5 obtained around the Hubble Deep Field-North (HDF-N), we made 
a composite spectrum of UV luminous (M1400 < —21.5 mag) LBGs at z ^ 5. The 
resultant spectrum shows a weak Lya emission and strong LIS absorptions which 
suggests that the bright LBGs at 2; ~ 5 have chemically evolved at least to ~0.1 solar 
metallicity. For a part of our sample in the HDF-N region, we obtained near-to-mid 
infrared data, which constraint stellar masses of these objects. With the stellar mass 
and the metallicity estimated from LIS absorptions, the metallicities of the LBGs at 
2; ~ 5 tend to be lower than those of the galaxies with the same stellar mass at 2^2, 



although the uncertainty is very large. 
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1. Introduction 

Spectroscopic studies of high-redshift (2; > 2 ~ 3) galaxies have been advanced recently. 
One of the largest spectroscopic samples of high-z galaxies is obtained by follow-up spectroscopy 
for Lyman Break galaxies (LBGs: e.g., Steidel et al. 2003) which are selected by using multi- 
band photometric data. So far, about a thousand of optical spectra have been obtained each for 
LBGs at z ~ 3 (e.g., Shapley et al. 2003) and star- forming galaxies at 2; ~ 2 (BM/BX objects: 
e.g., Steidel et al. 2004). Intensive spectroscopic studies based on magnitude- limited samples 
have also identified several hundreds of galaxies at z ~ 2 — 4 (e.g. FORS Deep Field (FDF) 
survey: Noll et al. 2004, VIRMOS-VLT Deep Survey: Le Fevre et al. 2005). 

These studies revealed spectroscopic features of galaxies at 2 ~ 2 — 3 in their rest-frame 
UV wavelength region, and deepened our understanding of their nature. Shapley et al. (2003) 
found relations between the strength of the Lya equivalent width (EW) and other spectroscopic 
features (e.g., EW of the low-ionized interstellar (LIS) metal absorption lines, UV continuum 
slope, velocity offset between Ly« and LIS absorptions) for the LBGs at z ~ 3, which sheds light 
on properties of an interstellar matter such as a gas outflow, a velocity dispersion of the gas, 
and a covering fraction of dust. With a part of FDF data, Mehlert et al. (2002) found that EW 
of the CIV AA1548, 1550 absorption line decreases from z ~ 1.4 to ~ 3.4. The decrease in the 
EW of CIV is confirmed with the full sample of FDF, and it seems to continue to z ~ 4 (Noll et 
al. 2004). Since the EW of CIV is an indicator of the metallicity in local starburst galaxies (e.g., 
Heckman et al. 1998), the result suggests a chemical evolution of galaxies in this redshift range 
(Mehlert et al. 2002), although the correlation between the metallicity and the EW of CIV 
has not been calibrated at high-redshift and thus is not ensured. Near-infrared (NIR) spectra 
of galaxies at 2 ~ 2 — 3 have also been obtained and give information about kinematics and 
metallicity of the galaxies. The velocity offset between nebula emission lines and interstellar 
absorption lines is generally seen in the LBGs, suggesting the presence of the gas outflow (e.g., 
Pettini et al. 2001). From the emissions lines, gas metallicity in the galaxies can be derived 
by using metallicity estimators (the i?23 and the N2 index); for example, Teplitz et al. (2000), 
Kobulnicky & Koo (2000), and Pettini et al. (2001) derived metallicity of LBGs at 2; ~ 3 which 
ranges typically from 0.1—0.9 solar value. Shapley et al. (2004) showed that majority of the 
seven NIR bright [Ksyega < 20 mag) BX galaxies have near the solar metallicity, and Erb et al. 
(2006) investigated the extended sample mainly with .ft's^yega < 21.5 mag and found the median 

* Based 011 data collected at Subaru Telescope, which is operated by the National Astronomical Observatory 
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metallicity of 87 galaxies at 2; ~ 2 to be about 0.5—0.6 solar value. 

Recently, surveys of galaxies at 2; ^ 5 have been carried out (e.g. Iwata et al. 2003; 
Lehnert & Bremer 2003; Stanway et al. 2003; Dickinson et al. 2004; Ouclii et al. 2004; Shimasaku 
et al. 2005; Yan et al. 2005; Bouwens et al. 2006; Yoshida et al. 2006), and challenges of finding 
galaxies even at 2; ~ 7 — 10 have also started (e.g., Pello et al. 2004; Kneib et al. 2004; Bouwens 
et al. 2005; Bouwens & Illingworth 2006; Mannucci et al. 2007). However, the progress of 
optical follow-up spectroscopy for galaxies at 2 ^ 4 — 5 is much slower than that for galaxies 
at 2 = 2^3 because objects are getting fainter, and characteristic features of UV spectrum of 
star- forming galaxies (e.g., Lya, LIS absorptions) redshift to the wavelength region where night 
sky emissions are very strong. Although attempts of spectroscopic observations of galaxies at 
2; ^ 5 were made (e.g., Spinrad et al. 1998; Dey et al. 1998; Weymann et al. 1998; Lehnert 
& Bremer 2003; Stanway et al. 2004), most of the galaxies were identified only with a Lya 
emission because of the low signal-to-noise ratio of their continuum light. 

We also started an optical spectroscopy of LBGs at 2 ~ 5. Targets were selected from 
our surveys for LBGs at 2 ~ 5 (Iwata et al. 2003, 2007) which consist of the deep and wide 
V,Ic, and z'-band imaging survey with the Subaru/ Suprime- Cam in the two target fields: the 
region including the Hubble Deep Field-North (the HDF-N region) and the J0053-I-1234 region. 
Our survey fields are quite wide (1290 arcmin^ in total), and we obtain 228 bright objects with 
z'j<^^ < 25.0 mag (L ^ L* in UV luminosity function (UVLF) of LBGs at z ^ 5). We reported 
an initial result of the spectroscopy for a part of the objects in the HDF-N region (Ando et al. 
2004: hereafter Paper I). Thanks to the object brightness and good observing conditions, we 
detected continuum features of the LBGs at 2 ~ 5 and identified the galaxies to be at 2; ~ 5 
without Lya emission line. The most striking result was that there is no or a weak (EWrest < 
lOA) Lya emission line in 7 bright LBGs identified. Inspired by the result, the deficiency of 
the strong Lya emission in the bright LBGs at 2 ~ 5 was investigated by using combined data 
of Paper I, our new results obtained in the J0053-I-1234 region, and those from the literature; 
the deficiency is generally seen in the LBGs at 2 ~ 5 and even at 2 ~ 6 (Ando et al. 2006). If 
the quenched Lya emission is due to the dust, this deficiency suggests that bright LBGs have 
larger amount of dust and higher metallicity than faint ones. This hypothesis suggests that 
bright LBGs experienced star formation earlier than that of faint ones: differential evolution 
of the high-redshift galaxies with their luminosity, although there are other possible origins of 
the deficiency (e.g., a velocity structure in/around the galaxy). 

In this paper, we present the data of spectroscopic observation of LBGs at 2 ~ 5 in the 
J0053-I-1234 region which were used in our previous paper (Ando et al. 2006). In Section 2, we 
describe a sample selection, observation, and data reduction. The features of obtained spectra 
of the LBGs at 2 = 4 — 5 in the J0053-I-1234 region are shown in Section 3. Combining this 
result with our previous results (Paper I), we described spectroscopic properties of bright LBGs 
at 2 ~ 5 and metallicity against stellar mass in the section. Throughout this paper, we adopt 
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flat A cosmology, Qm = 0.3, = 0.7, and Hq = 70 km s ^Mpc ^. The magnitude system is 
based on AB magnitude (Oke & Gunn 1983). 

2. Observation and Data Reduction 

The photometric sample of LBGs at z ~ 5 in the J0053+1234 region was selected based 
on a deep and wide broad-band (y,Ic, and z') imaging survey by using Subaru/Suprime-Cam 
(lye et al. 2004; Miyazaki et al. 2002): a total of 114 (236) candidates with z' < 25.0 (25.5) 
mag in an effective survey area (~800 arcmin^). Details of imaging observations and the color 
selection are described in Iwata et al. (2007). In order to detect continuum and absorption 
feature, we selected the bright {z' < 25.0) LEG candidates at z ~ 5 as the main spectroscopic 
targets. Because the entire survey field was too wide to obtain spectra for the whole sample 
of these LEG candidates with a limited field of view (6'0) of the Faint Object Camera and 
Spectrograph (FOCAS: Kashikawa et al. 2002) on the Subaru telescope and limited observing 
time, we selected a target MOS field to contain the bright LEGs as many as possible. When a 
slit of a bright target overlapped with a slit of another bright one, we chose the object with a 
higher central surface brightness. Finally, we designed one FOCAS mask which covers 5 bright 
LEGs. The mask also covered 8 fainter {z' > 25.0 mag) LEGs which are in the gaps of the slits 
of main targets on the mask. 

Spectroscopic observations with the Subaru/FOCAS were made on 2004 September 14 
under a clear condition, and a seeing during the observing run was typically ~ 0."8 (FWHM). 
We used a grism of 300 lines mm~^ blazed at 7500 A and the S058 order-cut filter (the same 
setting as Paper I), which gave wavelength coverage from 5800A to lOOOOA (depending on a 
slit position on the mask) with a pixel scale of 1.34A. The slit lengths were typically lO", and 
the slit widths were fixed to be 0."8, giving a spectral resolution of R~700 which was measured 
with night sky emission lines. An exposure time of each frame was 0.5 hours, and a total 
effective exposure time was 6 hours. 

The data were reduced by using IRAF^ with the same manner as Paper I except for the 
re-subtraction of the night sky emissions. We subtracted bias from object frames by using an 
overscan region of each frame and averaged bias frame in order to remove the characteristic 
bias-pattern of FOCAS CCDs. Flat-fielding was made with normalized average of dome-flat 
images. The wavelength calibration was made by using night-sky emission lines. The RMS 
error of the wavelength calibration was ~ 0.7A. After the sky subtraction, the spectrum of each 
object was carefully aligned and averaged. Eecause some objects show too low continuum S/N 
to subtract sky emission in each frame, we carried out sky subtraction after combining spectra 
for these objects. Five pixels were binned toward wavelength direction in order to improve 

^ Image Reduction and Analysis Facility, distributed by National Optical Astronomical Observatories, which 
arc operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 
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S/N, and one-dimensional spectra were obtained with APALL task of IRAF. Since our objects 
were faint, we did not fix the aperture width for the extraction but we determined it by eye 
for each object to trace it well. The flux calibration and the sensitivity correction were made 
for the final spectra by using standard stars of Wolfl346 and PG0823+546. The correction 
for the atmospheric A and B bands was also made by tracing those absorption features in the 
spectra of standard stars. For some of the objects, the sky subtraction was not well achieved 
due to the low-S/N of the continuum or the presence of the bad column of the CCDs. Thus 
we re-subtracted the scaled night sky emissions from the one- dimensional spectrum of such 
objects. 

3. Results and Discussions 

3.1. Redshift Determination 

Among 5 bright targets, we identified 2 objects to be LBGs at z = 4 — 5, and the resultant 
spectra are shown in Figure 1 (No.l: top left. No. 2: top right). LEG identifications were mainly 
made by the presence of the continuum depression shortward of the redshifted Lya due to the 
inter-galactic HI absorption and the Lya emission line near the Lya depression. In addition, 
the Lya emission line of each object shows an asymmetric shape which is considered to be 
characteristic feature of LBGs and Lya emitters (LAEs) at high- redshift. Note that emission- 
like features in the spectrum of object No.l at ~6300A and in the wavelength region longer 
than 7300A are the residuals of the night sky emissions since the sky subtraction was not 
well achieved due to the bad column in the CCD. For object No. 2, we identified some LIS 
absorption lines (Sill A1260, Ol+Sill A1303, and CII A1335) at almost the same redshift of 
the Lya which are characteristic features of nearby star- forming galaxies (e.g., Heckman et al. 
1998) and LBGs at z~3-5 (e.g., Steidel et al. 1996; Shapley et al. 2003; Frye et al. 2002; Paper 
I; Iwata, Inoue & Burgarella 2005), supporting the identifications of high redshift star-forming 
galaxies. Redshifts of these two bright LBGs were determined from the peak of Lya emission: 
4.797 and 4.267 for object No.l and No. 2, respectively. The error of the redshift determination 
was ~ 0.005. The rest of the bright sample were not identified to be at 2; ~ 5 because of their 
low S/N of the spectra. These objects tend to have lower central surface brightness compared 
with the confirmed objects. 

In addition to the bright targets, we identified 2 objects among faint bonus targets with 
z' > 25.0 mag to be LBGs at z = A — 5. Spectra of them are shown in Figure 1 (No. 3: bottom 
left, and No. 4: bottom right). These objects show a strong (EWobs ~ 220A and ~420A for 
object No. 3 and No. 4, respectively) single emission line; peak of the emission corresponds to 
the redshift of 4.391 and 4.491 for object No. 3 and No. 4, respectively if this emission is Lya. 
A clear asymmetry of the emission line profile and possible continuum depression shortward 
of the emission are seen in the spectrum of object No. 4, which supports this object is really 
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at z = 4 — 5. Because the redshift determinations for these objects were based only on the 
single emission line, we examine the possibility that they are foreground emission galaxies. 
For example, if object No. 3 (No. 4) is an [Oil] A3727 emitter at z = 0.76 (0.79), they should 
show the redshifted H/3 emission line at ~8550(8710)A and the [OIII]A5007 emission line at 
~8810(8970)A with EWobs of ~70(130)A and EWobs of ~110(210)A, respectively by assuming 
the median of the observed line ratios of star-forming galaxies at z = 0.3 ~ 1.0 (0.3 and 0.5, 
respectively: Kobulnicky & Kewley 2004). As shown in Figure 1, emission-like features are 
seen around at a part of the expected wavelengths in both spectra. However, these features are 
residual night-sky emissions, and there is no other significant (EWobs ^ lOOA) emission lines at 
the expected wavelengths. Thus these strong emissions are considered to be the Lya emission 
from LBGs at the target redshift. We can not completely rule out the possibility that they 
are foreground emission galaxies by only this analysis, because some of the expected lines of 
foregrounds are in the wavelength region of severe night sky emissions, and the line ratios have 
wide object-to-object variances (Kobulnicky & Kewley 2004). But the V — Iq colors of these 
objects (1.98 mag and 1.85 mag) are too red to regard them as star-forming galaxies at 2 < 1; 
we calculated the SEDs of star-forming galaxies at 2; < 1 by using the population synthesis code 
PEGASE version 2 (Fioc & Rocca-Volmerange 1997) with Salpeter IMF (Salpeter 1955) and 
by assuming constant star formation with an age of 100 Myr, which gives their V — Iq colors in 
observed-frame of —0.4 — 0.1 mag and 0.2 — 0.8 mag for E{B — V) = 0.0 mag and E{B — V) = 0.4 
mag, respectively with the dust extinction curve by Calzetti et al. (2000). Thus we consider 
both of the objects are at z = 4 — 5. 

3.2. Color Selection and Redshift Distribution of LBGs at z ^ 5 

Figure 2 shows the positions of observed objects in the V — Ic and Ic — z' two-color 
diagram from the present results and Paper I. Photometric properties of spectroscopically 
identified LBGs at 2: = 4 — 5 in the J0053+1234 region are summarized in Table 1. Filled circles 
show LBGs at 2; = 4 — 5 in this study: 2 bright objects as large circles and 2 faint objects as 
small ones. Filled squares show bright LBGs at 2 ~ 5 in our previous results (Paper I). Open 
triangles and open-inverse triangles show bright LBG candidates unidentified in our present 
observation and the previous results, respectively. So far, no bright spectroscopic targets (total 
of 22 LBG candidates) have been identified as foreground galaxies, though a half of them were 
not identified due to the low-S/N of their spectra. In order to examine our color selection 
criteria further, we also show the positions of the Galactic M stars which were identified in 
our previous observations as filled pentagons in Figure 2. Our selection criteria for LBGs at 
2 ~ 5 are reasonable to exclude interlopers, though more large and systematic spectroscopic 
observations are desirable. 

The redshift distribution of the spectroscopically confirmed LBGs at 2 ~ 5 in this study 
and Paper I is shown in Figure 3. In Figure 3, we plot the expected redshift distribution of our 
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Fig. 1. Spectra of LBGs at z = 4 - 5 in the J0053+1234 region. Top: Object No.l (left) and No.2 (right). 
Bottom: Object No.3 (left) and No. 4 (right). These spectra are smoothed with a boxcar over 3-pixeL 
Positions of hne features seen in LBGs at z = 3 and local starburst galaxies (e.g., Shapley et al. 2003; 
Heckman et al. 1998) are shown with vertical dashed lines. A scaled sky spectrum is shown in a lower part 
of each panel, and the atmospheric A-band and B-band absorptions are shown as vertical hatched regions. 

Table 1. Photometric properties of tiie LBGs identified to be at z ~ 5 in tlie J0053+1234 region. 



No. (ID) 


RA(J2000) 


DEC(J2000) 


z'* 


V-Ic^ 


Ic-z 


1 (106426) 


00:52:21.34 


+12:32:35.3 


24.07 


2.59 


0.23 


2 (104115) 


00:52:43.27 


+12:32:08.2 


24.22 


1.75 


0.16 


3 (091813) 


00:52:39.88 


+12:29:44.1 


25.03 


1.98 


0.06 


4 (093014) 


00:52:37.37 


+12:29:58.7 


25.26 


>1.85 


0.16 



* z' magnitude. MAG_AUTO from SExtractor (Bertin and Arnouts 1996) is adopted, 
t These values are measured with a 1".6 aperture. 
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Fig. 2. Positions of identified objects in the Ic — z' and V — Ic two-color diagram. Our color selection 
criteria for LBGs at 2 ~ 5 (Iwata et al. 2003; 2007) are indicated by thick lines. The objects confirmed to 
be at 2: ~ 5 in this study are shown as filled circles with each identification number in Table 1. The large 
circles show 2 bright objects, and small ones show 2 faint objects. Filled squares represent 7 bright LBGs 
at z 5 from Paper I. Open triangles and open-inverse triangles show bright LEG candidates unidentified 
in our present observation and the previous results (Paper I), respectively. Filled pentagons show objects 
identified to be foreground objects (Galactic M stars) from Paper I. A dashed (dot-dashed) line represents 
a color track of a model LEG spectrum with E{B — F) = 0.0 mag {E{B — V)= 0.4 mag) from Iwata et al. 
(2007). Small crosses are plotted on the tracks for z > 4.0 with an interval of 0.1. Model SEDs of LEGs 
are generated by using the population synthesis code PEGASE version 2 (Fioc & Rocca-Volmcrangc 1997) 
with Salpeter IMF (Salpeter 1955) and by assuming constant star formation with an age of 100 Myr. IGM 
absorptions at each redshift are calculated by using the analytic formula by Inoue et al. (2005). Calzctti's 
dust extinction curve (Calzetti et al. 2000) is assumed to make model spectra with E{B — V) = 0.4 mag. 
A dotted line refers to a color track of an elliptical galaxy (Coleman et al. 1980). Small open squares are 
plotted on the track with a redshift interval of 0.5. Small open pentagons indicate the colors of AO - M9 
stars calculated based on the library by Pickles (1998). 
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sample of LBGs at 2; ~ 5 as a solid line normalized at z = 4.7. The distribution is calculated 
based on the expected detection rates of the LBGs against the apparent magnitude and the 
redshift for each survey field (figure 6 by Iwata et al. (2007)). However, the detection rates 
were derived through simulations in the two-color diagram to estimate the completeness of 
the photometric catalog, and they do not include the difference of survey volumes for each 
redshift bin and the difference of the number density at the corresponding UV magnitude for 
an apparent magnitude at each redshift bin. In order to derive the expected distribution and 
compare it with the observed redshift distribution of the spectroscopic sample, we corrected 
the detection rates by adopting the redshift bin size of Az = 0.2 and assuming the UVLF of the 
z ~ 5 LBGs (Iwata et al. 2007). Since almost all our spectroscopic sample covers z' magnitude 
brighter than 25.5 mag, we calculated the expected redshift distributions for three magnitude 
ranges of 24.0—24.5 mag, 24.5—25.0 mag, and 25.0—25.5 mag for each field, and took the 
average of them. The distribution of spectroscopically confirmed LBGs is broadly consistent 
with the expected distribution, but there seems to be a deficiency of the LBGs at around ^ ~5.0. 
Although the sample size is still small and the deficiency may not be significant, this may be 
due to the presence of severe night sky emission lines which prevent us from detecting the 
continuum depression and/or Lya emission. The night sky emission lines whose wavelengths 
are converted to the redshifts corresponding to the redshift of Lya are also shown with dotted 
line in the figure; the sky emission lines are severe for the redshifts larger than 4.9. 

As described in section 3.1, spectroscopically identified LBGs in the J0053 region have 
redshifts of 4.3 — 4.8 with a mean of 4.5, which is slightly lower than that of the objects in Paper 
I with redshifts of 4.5 — 5.2 (mean value of 4.7). In fact, as seen in Figure 2, the 4 objects in the 
J0053 region are bluer in the V — Ic color than the 7 objects in Paper I. The mean V — Ic color 
of spectroscopic sample is 2.62 mag and 2.04 mag for the objects in the HDF-N region and in 
the J0053 region, respectively. However the difference may be due to the variation among small 
numbers of objects, and we should have more objects with spectroscopic redshifts to examine 
this point further. 

3.3. Spectra of LBGs at z 5 in the J0053+1234 region 

As shown in Figure 1, spectra of bright LBGs in the J0053+1234 region show Lyct and 
LIS absorptions (Sill A1260, Ol+Sill A1303, and CII A1335) which are typical line features 
of star-forming galaxies at low and high redshifts. The Lya EWrestS of object No.l and No. 2 
are 14±3A and 21±4A, respectively, and their observed fluxes are 2.5 xlO^^'^ erg s~^ cm~^ and 
2.4x10"^'' erg s~^ cm~^, respectively corresponding to the luminosities of 5.8x10^^ erg s~^ and 
4.3x10^^ erg s~^, respectively. These values of Lya do not include absorption component and 
are not corrected for the IGM absorption. The errors of EW values in this paper include the 
uncertainty in determination of the continuum level and the RMS error of the continuum. The 
values of EW of Lya are larger than the average value of our previous results for 7 bright LBGs 
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Fig. 3. Redshift distribution of LBGs at z ~ 5 from this study and our previous spectroscopic results 
(Paper I). Tlie night sky emission lines whose wavelengths are converted to the redshifts corresponding to 
the redshift of Lya are also shown with dotted line. Expected redshift distribution normalized at z = 4.7 
for LBGs at z ^ 5 is also plotted as solid line (right vertical axis, see text for the details). 

(2.5±3.7A, Paper I), but they are not so large to be detected as definitive LAEs. Typical Lya 
EWrest value for selecting LAEs is 15~20A. In contrast to the result for the bright targets, the 
Lya emissions of the faint LBGs are strong; EWrost of object No. 3 and No. 4 are 40±12A and 
77±23A, respectively. The observed fluxes of Lya are 3.7x10"^^ erg s~^ cm~^ and 8.3x10"^^ 
erg s~^ cm~^, respectively corresponding to the luminosities of 7.1x10^^ erg s~^ and 1.7x10^^ 
erg s~^, respectively. We summarize the observed quantities in Table 2. 

The bright LBGs at z ~ 5 tend to show no or weak Lya emission; average EWrest of 
Lya emission of the total of 9 bright LBGs is ~ 5.8 ± 7.2A, while strong Lya emission is seen 
among the faint LBGs. Because we can not deny the presence of faint LBGs with no or weak 
Lya emission from the current data, these results suggest, at least, a deficiency of strong Lya 
emission in bright LBGs at 2; ~ 5 and might suggest the UV luminosity dependence of the Lya 
emission. Combining with other spectroscopic results from the literature, we found that the 
deficiency is generally seen in LBGs at z ~ 5 and possibly in LBGs at z ~ 6 (Ando et al. 2006). 
The trend is also seen in the galaxies at z = 4^6 (Vanzella et al. 2006; 2007). Possible origins 
of the deficiency of strong Lya emission are discussed by Ando et al. (2006). The deficiency or 
the UV luminosity dependence of the Lya are also seen in the recent results of spectroscopy of 
LBGs at z ~ 3; Shapley et al. (2006) showed that UV luminous [Muy ~ —21.6 mag) LBGs at 
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Table 2. Spectroscopic properties of the LBGs identified to be at 2 = 4 — 5 in the J0053+1234 region. 



No. (ID) 


Redshift* 


EW(Lya)t 


Ft * 






1 (106426) 


4.797 


14±3 


2.5x10- 


17 


5.8x10^^2 


2 (104115) 


4.267 


21±4 


2.4x10" 


17 


4.3xl0'*2 


3 (091813) 


4.391 


40±12 


3.7x10- 


17 


7.1x10^2 


4 (093014) 


4.491 


77±23 


8.3x10- 


17 


1.7x10^3 



' Kedshitts are determined from the peak ot Lya emission. I'he error ot the redshift 
determination is ~ 0.005. 

t Rest-frame equivalent width of Lyo emission in units of angstroms, not including 
Lyo absorption. EWs are taken to be positive for emission lines. The error is 
estimated from the uncertainty of the continuum level determination and the RMS 
error of the continuum. 

t Observed flux of the Lya emission line in units of erg s~^ cm~^. 

§ Luminosity of the Lya emission line in units of erg s~^. The slit correction is not 

applied to the values. 

z ~ 3 tend to show weaker Lya emission than the total LBG sample of Shapley et al. (2003). 
From the re-binned sample of Shapley et al. (2003) with their UV luminosity, Keel (2006) 
showed that Lya emission is weaker in the UV luminous {Muv < —21.0 mag) LBGs than in 
fainter LBGs. 

For object No. 2, we identified LIS absorption lines (Sill A1260, 01+Sill A1303, and 
Gil A1335) at almost the same redshift as the Lya emission. We also identified SilV A1394 
(EWrest = -2.0 ±1.1 A) 2 and GIV AA1548,1551 (EW^cst = -4.9±1.4A) which are also seen 
in the spectra of LBGs at ~ 3 (e.g., Shapley et al. 2003). The average EW of the three LIS 
absorption lines is — 4.2 ± 1.4A which is larger than the average value of our previous result 
(— 2.8± I.2A : Paper 1). The peak of Lya emission line is redshifted by 280±210 km s~^ relative 
to the average of the LIS lines. The peak-to-valley velocity offset which may be related to a 
large scale gas outflow was also reported in LBGs at 2; = 3 ~ 5 (e.g., Shapley et al. 2003; Frye 
et al. 2002; Paper I). The offset value of object No. 2 is smaller than typical value of LBGs at 
z ~ 3 (500 - 600 km s"^: Shapley et al. 2003) and that of LBGs at 2 ~ 5 (~ 600 km s"^: Paper 
I). For object No.l, we could not identify significant absorption features due to the residuals of 
the sky subtraction. 

3.4- Composite Spectrum of Luminous LBGs at z ^ 5 

In order to overcome the low-S/N of each spectrum and investigate typical spectroscopic 
properties of luminous LBGs at z ~ 5, we produced a composite spectrum of 8 LBGs at z ~ 5 
with Mi4oo < —21.5 mag by averaging the spectra obtained in this study and our previous 

2 We resolved SilV AA1394,1403 doublet, but A1403 component was not detected significantly (comparable to 
Icr), which is consistent with the observed ratio for LBGs at z 3 (Shapley et al. 2003). 
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study. We derived the absolute UV magnitude from the observed broad band magnitude to the 
rest-frame 1400A magnitude assuming a continuum slope /? {fx oc A'^) of —1 which is a typical 
value for LBGs at z ~ 3 (Shapley et al. 2003). Before co-adding the spectra, all wavelengths 
were converted to the rest-frame scale using individual redshifts, and each flux density was 
normalized at 1250A. 

Figure 4 shows the resultant composite spectrum re-binned to common dispersion of lA 
pixel"^. Note that the emission-like features except for the Lya (e.g., around 1275A, 1425A) 
come from the residuals of sky emission lines, and these features are not real. The continuum 
depression shortward of Lya by IGM is clearly seen in the composite spectrum. A depression 
factor Da (Oke & Korycansky 1982) for the composite spectrum is ~0.6, which agrees with 
the value of the composite spectrum of LBGs at 2; ~ 5 (Paper I) and the values of QSOs at 
2; ~ 5 (0.1 — 0.7: e.g., Songaila & Cowie 2002; Songaila 2004). As expected from the previous 
subsection, EWrest of the Lya emission is not large (7.6±0.8A), and LIS absorptions are also 
seen clearly. The EWs of LIS absorptions are -2.7±0.36A, -1.8±0.30A, and -1.5±0.29A 
for Sill A1260, Ol+Sill A1303, and CII A1335, respectively, which suggests that bright LBGs 
have been chemically polluted to some extent by metal at z ~ 5. The SilV AA1394,1403 is also 
seen in the composite spectrum, and the EWrestS of the SilV A1394 and the SilV A1403 are 
-2.0 ±0.32 A and -1.0 ±0.26 A, respectively. 

The EWrest of Lya is smaller than that of the composite spectrum of LBGs at 2; ~ 3 
(15. lA: Shapley et al. 2003); the EW values of LIS absorptions are comparable to or somewhat 
larger than those of z ~ 3 (-1.7A, -2.3A, and -I.SA for Sill A1260, OI±SiII A1303, and CII 
A1335, respectively). However, the sample luminosity range of our LBGs at 2 ~ 5 {z' ^ 24.8 
mag corresponding to M1400 ^ —21.5 mag) is different from that of LBGs at 2: ~ 3 (mainly 
i? = 23.5 — 25.5 mag corresponding to M1400 = —21.8 ~ —19.8 mag). As described in the previous 
subsection, we found the deficiency of the strong Lya emission in luminous LBGs at 2 ~ 5 
and pointed out the possible presence of the UV luminosity dependence on the spectroscopic 
features. Thus we can not discuss the evolution of the spectroscopic feature of LBGs from a 
simple comparison of the composite spectra. 

Although the low S/N around the wavelength range longer than 1400A does not allow us 
to estimate metallicity by using 1425 index (e.g., Rix et al. 2003) and CIV index (e.g., Mehlert 
et al. 2002), we roughly estimate mean metallicity of the 8 luminous LBGs at 2 ~ 5 from LIS 
absorption lines of the composite spectrum assuming the relation between the mean EW of the 
three LIS absorptions (Sill A1260, OI±SiII A1303, and CII A1335) and a metallicity calibrated 
with local star-forming galaxies (Heckman et al. 1998). Since the sample by Heckman et al. 
(1998) does not cover the low metallicity (i.e., weak LIS absorptions) region, we extrapolated 
the EW(LIS)-metallicity relation in the middle-to-high (12±log(O/H)^8.0) metallicity region 
to the lower metallicity region. The estimated metallicity of 12±log(0/H) is about 7.6±0.6 
corresponding to ~0.1 (0.3—0.02) solar value (8.66: Asplund et al. 2004). The errors come 
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Fig. 4. Composite spectrum of 8 luminous (M1400 < —21.5 mag) LBGs at z 5 from this study and 
Paper I. Note tliat most of emission-like features (e.g. around 1275A, 1425A) come from the residuals of 
sky emission lines. 

from the uncertainty of the EW measurement as well as that of the empirical relation between 
the EW(LIS) and metallicity. It should be worth noting that the LIS lines are saturated, and 
EW(LIS) reflects not only the metallicity but also the velocity structure and covering fraction 
of the gas (e.g., Heckman et al. 1998; Shapley et al. 2003; Noll et al. 2004). Erb et al. (2006) 
show a composite rest-frame UV spectrum of objects at z ~ 2 with a mean stellar mass of 
7 X IO^^Mq and mean metallicity of 12+log(0/H)~ 8.55 derived from N2 index method. The 
EW(LIS) in the spectrum is ~ 2 — 3A corresponding to 12-|-log(0/H) of ~ 7.6 — 8.2 if we use the 
local relation by Heckman et al. (1998). The result suggests that the local EW(LIS)-metallicity 
relation may systematically underestimate the metallicity of galaxies at high redshift. 

3.5. A Hint for Chemical Evolution of massive LBGs at z ^ 5? 

We derived the individual metallicity for 8 LBGs at 2; ~ 5 with the LIS absorption in 
this study and Paper I with the same manner mentioned above. For the 3 objects in the HDF- 
N region among the 8 LBGs, Subaru/CISCO K' band data and Spitzer/IRAC mid-IR data 
are available, and we can estimate their stellar masses by fitting their rest-frame UV to optical 
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SEDs to the model SEDs. (Details of NIR and MIR data of our LBGs at 2 ~ 5 will be described 
elsewhere.) The metallicities (12+log(0/H)) of 104268 (No.l of Paper I), 144200 (No.2 of Paper 
I), and 148198 (N0.8 of Paper I) are 7.7±0.8, 7.7±0.8, and 8.7±0.9, respectively. For the stellar 
mass estimation, we adopted the same manner used in Iwata, Inoue & Burgarella (2005). We 
used the PEGASE version 2 population synthesis model (Fioc & Rocca-Volmerange 1997) with 
Salpeter IMF (Salpeter 1955), constant star formation history, an the Calzetti's dust extinction 
curve (Calzetti et al. 2000) to make model spectra. IGM absorptions were calculated for the 
given redshift by using the analytic formula by Inoue et al. (2005). Then the observed SEDs 
were fit to model SEDs, and the best fit model was obtained by minimizing . Resultant 
stellar masses (best fit value with 68% confidence level) of the objects 104268, 144200, and 
148198 are 5.82^1^ x lO^^M©, 3.33l^|t x lO^^M©, and LSS^olo x lO^M©, respectively. 

In Figure 5, open circles represent the stellar masses and the metallicities of the three 
luminous LBGs at 2;~5, and filled circle represents the average value of them. Pluses, diamonds, 
and triangles show the data for galaxies at 2~0.1, 2;~0.7 and 2; ~ 2, respectively (Tremonti 
et al. 2004; Savaglio et al. 2005; Erb et al. 2006), which show the presence of mass-met allicity 
relation at each redshift and its evolution from 2~2to2;~0^. Although the uncertainties 
are quite large (and may contain some systematics), and our way of metallicity estimation is 
different from those for galaxies at other redshifts, the metallicities of the LBGs at 2; ~ 5 tend 
to be lower than those of the galaxies with the same stellar mass at z ^ 2 ^. 

We are grateful to the FOCAS team, especially Youichi Ohyama, and the staffs of 
Subaru telescope for their supports during our observation. We also thank Kiyoto Yabe and 
Marcin Sawicki for useful discussions in deriving stellar masses. We appreciate the anonymous 
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Center (ADAC) and Subaru Telescope of the National Astronomical Observatory of Japan. 
MA is supported by a Research Fellowship of the Japan Society for the Promotion of Science 
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Society for the Promotion of Science (17540216). This work is supported by the Grant-in-Aid 
for the 21st Century COE "Center for Diversity and Universality in Physics" from the MEXT 
of Japan. 



The metallicity of galaxies at z 2 (Erb ct al. 2006) was derived from N2 index (e.g., Pettini & Pagel 2004) 
which is different from methods used by Tremonti et al. (2004) and Savaglio et al. (2005) and may have a 
problems in calibration (e.g., Erb et al. 2006) which may partly be a cause for the offset. 
The stellar masses by Erb et al. (2006) were derived by using Chabrier IMF. Thus our data points should 
be ~1.8 times smaller (e.g., Erb et al. 2006) in Figure 5 in order to compare them with the data by Erb et 
al. (2006). 
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Fig. 5. Stellar mass-metallicity relation for galaxies at various redshifts. Open circles represent our results 
of 3 luminous LBGs at 2; ~ 5, and filled circle represents the average value of them. Pluses, diamonds, 
and triangles show the results for galaxies at z ~ 0.1, z ~ 0.7 and z ~ 2, respectively (Tremonti et al. 
2004; Savaglio et al. 2005; Erb et al. 2006). 
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